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We analyze the prospects of a feasible, Brookhaven National Laboratory based, very long baseline 
(BVLB) neutrino oscillation experiment consisting of a conventional horn produced low energy wide 
band beam and a detector of 500 kT fiducial mass with modest requirements on event recognition 
and resolution. Such an experiment is intended primarily to determine CP violating effects in the 
neutrino sector for 3-generation mixing. We analyze the sensitivity of such an experiment. We 
conclude that this experiment will allow determination of the CP phase Sep and the currently 
unknown mixing parameter 813, if sin 2 2#i3 > 0.01, a value ~ 15 times lower than the present 
experimental upper limit. In addition to #13 and Sep, the experiment has great potential for precise 
measurements of most other parameters in the neutrino mixing matrix including Am^ , sin 2 2623 , 
Atria x sin 2#i2 , and the mass ordering of neutrinos through the observation of the matter effect in 
the Vp, — > appearance channel. 
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Introduction 

Measurements of solar and atmospheric neutrinos have 
provided strong evidence for non-zero neutrino masses 
and mixing [1-3]. Atmospheric results have been fur- 
ther strengthened by the K2K collaboration's accelera- 
tor based results [4]. The Solar neutrino results have 
been confirmed by the KamLAND collaboration in a re- 
actor based experiment that has shown that the large 
mixing angle (LMA) solution is most likely the correct 
one [5]. Interpretation of the experimental results is 
based on oscillations of one neutrino flavor state, v^v^ 
or v T , into the others, and described quantum mechani- 
cally in terms of neutrino mass eigenstates, V\, and 1^3. 
The mass squared differences involved in the transitions 
are measured to be approximately Ato^i = m { v 2) 2 — 



m(i'i) 2 = (5 — 20) x 10 5 cV 2 for the solar neutrinos and 
Am| 2 ee m(v 3 f - m(y 2 f = ±(1.6 - 4.0) x 10~ 3 cV 2 for 
the atmospheric neutrinos, with large mixing strengths, 
sin 2 2#i2 « 0.86 and sin 2 26> 2 3 ~ 1.0 in both cases. 

These existing data on neutrino oscillations [1-3] and 
the prospect of searching for CP violation make clear 
that the next generation of oscillation experiments must 
be significantly more ambitious than before. In particu- 
lar, the source of neutrinos needs to be accelerator based 
so that both the neutrino flavor content and the energy 
spectrum of the initial state can be selected. Several 
approaches have been explored in the literature. These 
involve either a narrow band "off-axis" beam produced 
with a conventional magnetic focusing system [6-8] or a 
neutrino factory based on a muon storage ring [9]. In 
this paper, we show that the currently favored param- 
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eters open the possibility for an accelerator based very 
long baseline (BVLB) experiment that can explore both 
solar and atmospheric oscillation parameters in a single 
experiment, complete the measurement of the mixing pa- 
rameters, and search for new physics. 

The Experimental Strategy 

There are four measurements of primary interest that 
can be addressed with the experimental arrangement de- 
scribed in this paper. Using the parameter convention 
described in [1], they are: 

(i) Definitive observation of oscillatory behavior 
(multiple oscillations) in the disappearance 
mode and precise determination of the oscillation 
parameters Am§ 2 and sin 2 2£> 23 , with statistical er- 
rors ~ ±1%. 

(ii) Detection of the oscillation — > v e in the ap- 
pearance mode and measurement of the parame- 
ter sin 2 2#i3. This will involve matter enhancement 
and also allow definitive measurement of the sign 
of Am§ 2 , i-e., which neutrino v% or is heavier. 

(iii) Measurement of the CP violation parameter, 
Sep, in the lepton sector. 

(iv) Measurement of Am^ x sin2#i 2 in the appear- 
ance channel of — ► v e . 

We describe how these measurements could be car- 
ried out with good precision in a single accelerator based 
experiment. For precise and definitive measurement of 
oscillations we must observe multiple oscillation nodes 
in the energy spectrum of reconstructed muon and elec- 
tron neutrino charged current events. The multiple node 
signature is also necessary in order to distinguish be- 
tween oscillations and other more exotic explanations 
such as neutrino decay [10] or extra dimensions [11] 
for the muon neutrino deficit in atmospheric neutrinos. 
Since the cross-section, Fermi motion, and nuclear effects 
limit the statistics and the energy resolution (for recon- 
structed neutrino energy) of low energy charged current 
events, we use neutrinos with energies greater than few 
hundred MeV and clean events with a single visible muon 
or electron that are dominated by quasielastic scattering, 
i>l + n — > I + p, where I denotes a lepton. Figure 1 shows 
that the distance needed to observe at least 3 nodes in the 
reconstructed neutrino energy spectrum is greater than 
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FIG. 1: (color) Nodes of neutrino oscillations for disappear- 
ance (not affected by matter effects) as a function of oscilla- 
tion length and energy for Am| 2 = 0.0025 eV 2 . The distances 
from FNAL to Soudan and from BNL to Homestake are shown 
by the vertical lines. 



2000 km for Am| 2 = 0.0025 eV 2 , the currently favored 
value from Super-Kamiokande atmospheric data [2]. A 
baseline of longer than 2000 km coupled with a wide band 
beam with high flux from 0.5 to 6 GeV will provide a 
nodal pattern in the — > disappearance channel and 
good sensitivity over a broad range of Am| . 

The low energy wide band beam and the very long 
baseline also have a number of important advantages for 
the appearance channel of — > v e . These advantages 
can be summarized using Figure 2 which shows the prob- 
ability of — > v e oscillation as a function of neutrino 
energy for the distance of 2540 km, the distance between 
BNL and the Homestake mine in South Dakota. The os- 
cillation parameters that we have assumed are indicated 
in the figure and the caption. We define the natural 
mass ordering (NO) of neutrinos to be 771 3 > m 2 > mi, 
and the other two possibilities (m 2 > mi > 7713 and 
mi > to 2 > m^) are collectively called reversed order- 
ing (RO). Both RO possibilities have very similar phys- 
ical consequences, but mi > ?n 2 > 7713 is strongly dis- 
favored by the explanation of the Solar neutrino deficit 
using matter enhanced neutrino oscillation (the LMA so- 
lution). Since neutrinos from an accelerator beam must 
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FIG. 2: (color) Probability of — » v £ oscillations at 2540 
km. The calculation includes the effects of matter. The dot- 
ted (Sep = 45°) and solid (Sap = 0°) curves are for NO and 
the lower dot-dashed (Sep = 0°) curves is for RO. The param- 
eters used for the figure are sin 2 2#i2 = 0.86, sin 2 2#23 = 1.0, 
and sin 2 26>i 3 = 0.04 and Am|j = 7.3 x 10~ 5 el/ 2 , Am 2 2 = 
0.0025 eV 2 . 



pass through the Earth to arrive at a detector located 
2540 km away, the probability in Figure 2 includes the ef- 
fects of matter, which enhance (suppress) the probability 
above 3.0 GeV for NO (RO) [12]. Therefore the appear- 
ance probability above 3.0 GeV is sensitive to both the 
mass ordering and the parameter sin 2 26*13 . The prob- 
ability in the region 1.0 to 3.0 GeV is less sensitive to 
matter but much more sensitive to the CP phase Sep as 
shown in Figure 2 [13]. The increase in the probability be- 
low 1.5 GeV is due to the presence of terms involving the 
solar oscillation parameters, Am^ and sin 9i%. There- 
fore, the spectrum of electron neutrino events measured 
with a wide band beam over 2540 km with sufficiently 
low background has the potential to determine sin 2 26*13, 
Sep, and the mass ordering of neutrinos as well as Am^ 
because these parameters affect different regions of the 
energy spectrum. In the following, we examine how well 
the parameters can be determined and the implications 
for the detector performance and background. 



Proposed Experimental Configuration 

The high energy proton accelerator, to be used for 
making the neutrino beam, must be intense (~ 1 MW in 
power) to provide sufficient neutrino-induced event rates 
in a massive detector very distant from it. Such a long 
baseline experimental arrangement [14] may be realized 
with a neutrino beam from the upgraded 28 GeV pro- 
ton beam from the Alternating Gradient Synchrotron at 
the Brookhaven National Laboratory (BNL) and a water 
Chercnkov detector with 0.5 megaton of fiducial mass 
at cither the Homcstake mine in South Dakota or the 
Waste Isolation Pilot Plant (WIPP) in New Mexico, at 
distances ~2540 km and ^2900 km from BNL, respec- 
tively A version of the detector is described at length 
in [15] and is a candidate for location in the Homcstake 
Mine, where it would occupy five independent cavities, 
each about 100 kton in fiducial volume. Another version 
of the detector in a single volume is described in [16]. For 
much of the discussion below we have used the distance 
of 2540 km to Homcstake for our calculations. Although 
the statistics obtained at 2900 km to WIPP will be some- 
what smaller, many of the physics effects will be larger, 
and therefore the resolution and the physics reach of the 
BVLB experiment are approximately independent of the 
choice between Homcstake or WIPP. 

The accelerator upgrade as well as the issues regard- 
ing the production target and horn system are described 
in [17]. We briefly describe it here for completeness. 
Currently, the BNL- AGS can accelerate ~ 7x 10 13 pro- 
tons upto 28 GeV approximately every 2 sec. This cor- 
responds to average beam power of about 0.16 MW. 
This average power could be upgraded by increasing 
the repetition rate of the AGS synchrotron to 2.5 Hz 
while keeping the number of protons per pulse approx- 
imately the same. Currently a 200 MeV room temper- 
ature LINAC in combination with a small synchrotron, 
called the Booster, injects protons into the AGS at 1.2 
GeV. The process of collecting sufficient number of pro- 
tons from the Booster into the AGS takes about 0.6 
sec. Therefore, for 2.5 Hz operation the Booster must 
be replaced by a new injector. A new super-conducting 
LINAC to replace the Booster could serve the role of a 
new injector; the remaining modifications to the AGS 
are well understood and they involve power supplies, the 
RF system, and other rate dependent systems to make 
the accelerator ramp up and down at 2.5 Hz. We ex- 
pect the final upgraded accelerator configuration to yield 
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8.9 x 10 13 protons every 400 ms at 28 GeV. Our stud- 
ies show that this intensity can be further upgraded to 
2 MW by incremental improvements to the LINAC and 
the AGS repetition rate. 

In addition to the 1 MW accelerator upgrade, the pion 
production target, the horn focusing system, and the de- 
cay tunnel, which must be aimed at an angle of 11.26 
degrees into the earth (towards Homestake) need to be 
built. Our current studies indicate that a carbon-carbon 
composite target inserted into a conventional aluminum 
horn cooled with forced Helium and water will function 
in the 1 MW beam for a sufficient length of time. Our 
current plan is to build a hill instead of an underground 
tunnel to accommodate the 200 m long pion decay tun- 
nel. We consider both the rapid cycling accelerator up- 
grade which maintains a relatively low intensity per pulse 
and the 1 MW capable target technically less risky than 
other alternatives. For comparison, at JPARC [18] the 
first (preliminary) stage of 0.75 MW is to be achieved at 
~ 1/3 Hz with 3.3 x 10 14 protons per pulse at 50 GeV. 
At Fermilab a 2 MW upgrade to the 120 GeV main in- 
jector calls for cither an 8 GeV super-conducting LINAC 
or a new rapid cycling proton synchrotron as injectors 
for the main injector[19]. As explained in the following, 
the combination of the very long baseline and a 500 kT 
detector allows the BVLB experiment to use a lower av- 
erage power at 1 MW compared to the proposed JPARC 
or NuMI(Off axis) [7] projects and achieve wider physics 
goals. This reduces the technical difficulties (liquid mer- 
cury or moving targets, radiation damage, shielding of 
personnel, etc.) associated with using beam powers in 
excess of 1 MW. 

We have performed detailed simulations of a wide band 
horn focused neutrino beam using 1 MW of 28 GeV pro- 
tons on a graphite target. The neutrino spectrum ob- 
tained in these simulations and used for the results in this 
paper is shown in Figure 3 normalized for POT (protons 
on target). We calculate that with such a beam (total 
intensity of ~ 4.7 x 10~ 5 v/m 2 /POT at distance of 1 km 
from the target) we will obtain about 60,000 charged cur- 
rent and 20,000 neutral current events in a total exposure 
lasting 5 x 10 7 seconds in a 0.5 megaton water Cherenkov 
detector located at the Homestake mine 2540 km away 
from BNL. For the results reported below, we use the 
particle spectra obtained by making cuts to select single 
muon or electron quasielastic events. Events with multi- 
ple particles (about 4 times higher in statistics) could be 
used to further enhance the statistical significance of the 
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FIG. 3: (color) The simulated wide band neutrino flux for 28 
GeV protons on a graphite target used for the calculations in 
this paper. (POT=protons on target). 



effects. 

A conventional horn focused beam can be run in either 
the neutrino or the anti-neutrino mode. In this paper we 
show that most of the physics program in the case of 3- 
generation mixing can be carried out by taking data in 
the neutrino mode alone. It is well known that exper- 
iments that must take data with anti-neutrinos, which 
have lower cross section, must run for long times to accu- 
mulate sufficient statistics [7] . This problem is largely al- 
leviated in our method. We will discuss the anti-neutrino 
beam in a separate paper; here we wish to focus on the 
physics reach of the neutrino beam alone [20]. 



f M Disappearance 

An advantage of the very long baseline is that the 
multiple node pattern is detectable over the entire al- 
lowed range of Am 2 2 . The sensitivity to Am^ ~ 
1.24 eV 2 extends to the low value of about 5 x 

L[kmJ 

10 _4 eV 2 , significantly below the allowed range from Su- 
per Kamiokandc. At energies lower than ~ 1 GeV, the 
energy resolution will be dominated by Fermi motion and 
nuclear effects, as illustrated in Figure 1. The contribu- 
tion to the resolution from water Cherenkov track recon- 
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FIG. 4: (color) Spectrum of expected single muon events in a 
0.5 MT water Cherenkov detector at 2540 km from the source. 
We have assumed 1 MW of beam power and 5 x 10 7 sec of 
data-taking. The top histogram is without oscillations; the 
middle histogram with error bars is with oscillations. Both 
histograms include the dominant single pion charged current 
background. The bottom histogram shows this background 
contribution to the oscillating spectrum. This plot is for 
Am§ 2 = 0.0025 eV 2 . The error bars correspond to the sta- 
tistical error expected in the bin. At low energies the Fermi 
movement, which is included in the simulation, will dominate 
the resolution. 

struction depends in first approximation on the photo- 
multiplier tube coverage. With coverage greater than 
10%, a reconstruction energy resolution of better than 
~ 10% should be achieved [6]. The simulated spec- 
trum of the expected disappearance signal includ- 
ing backgrounds and resolution is shown in Figure 4 
for Am§ 2 = 0.0025 eV 2 as a function of reconstructed 
neutrino energy. The background, which will be pri- 
marily charged current events, will also oscillate and 
broaden the dips in the nodal pattern. From this spec- 
trum we estimate that the determination of Am§ 2 will 
have a statistical uncertainty of approximately ±0.7% at 



Am| 2 = 0.0025 eV 2 and sin 2 26> 23 = 1. The experiment 
can determine sin 2 26* 2 3 > 0.99 at 90% confidence level. 
Within the parameter region of interest there will be little 
correlation in the determination of Am§ 2 and sin 2 2# 2 3. 
The precision of the experiment is compared in Figure 5 
with the precision expected from MINOS [21] and the re- 
sult from Super-Kamiokandc [2]. Since the statistics and 
the size of the expected signal (distortion of the spec- 
trum) are both large in the disappearance measurement, 
we expect that the final error on the parameters will be 
dominated by the systematic error. A great advantage 
of the very long baseline and multiple oscillation pat- 
tern in the spectrum is that the effect of systematic er- 
rors on flux normalization, background subtraction, and 
spectrum distortion due to nuclear effects or detector cal- 
ibration is small. The error on the overall detector energy 
scale is expected to be the dominant systematic error [14] . 
The large event rate in this experiment will allow us to 
measure Am| 2 precisely in a short period of time; this 
measurement will be important to predict the shape of 
the appearance signal which we now discuss. 



— > i/ e Appearance 

The importance of matter effects on long-baseline neu- 
trino oscillation experiments has been recognized for 
many years [12, 22]. For our study, we have included 
the effect of matter with a full numerical calculation tak- 
ing into account a realistic matter profile in the earth, 
following [23]. For our present discussion, it is useful to 
exhibit an approximate analytic formula for the oscilla- 
tion of — > v e for 3-generation mixing obtained with 
the simplifying assumption of constant matter density 
[24, 25]. Assuming a constant matter density, the oscil- 
lation of — > v e in the Earth for 3-generation mixing is 
described approximately by Equation 1. In this equation 
a = Am^/Am 2 .!, A = Am| 1 £/4£', A = 2V A £'/Am| 1 , 
V = \/2Gpn e - n e is the density of electrons in the Earth. 
Recall that Am^ = Am 2 2 + Ato 21 . Also notice that 
AA = LGpTiel 'y/2 is sensitive only to the sign of Ato 2 ^. 
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For anti-neutrinos, the second term in Equation 1 has the opposite sign. It is proportional to the following CP 
violating quantity. 



J cp — sin 6*i2 sin #23 sin 6*13 cos 6*12 cos #23 cos 2 6*13 sin 5c p 

I 



(2) 



Equation 1 is an expansion in powers of a. The ap- 
proximation becomes inaccurate for Aml 2 L/AE > it/2 
as well as a ~ 1. For the actual results shown in Figure 
2 we have used the exact numerical calculation, accurate 
to all orders. Nevertheless, the approximate formula is 
useful for understanding important features of the ap- 
pearance probability: 1) the first 3 terms in the equa- 
tion control the enhancement (for NO) or suppression 
(for RO) of the oscillation probability above 3 GeV; 2) 
the second and third terms control the sensitivity to CP 
in the intermediate 1 to 3 GeV range; and 3) the last 
term controls the sensitivity to Am^ at low energies. 

While the disappearance result will be affected 
mainly by systematic errors, the — > v e appearance 
result will be affected by the backgrounds. The v e signal 
will consist of clean, single electron events (single show- 
ering rings in a water Chcrcnkov detector) that result 
mostly from the quasielastic reaction v e + n — > e~ + p. 
The main backgrounds will be from the electron neutrino 
contamination in the beam and reactions that have a tt° 
in the final state. The 7r° background will depend on how 
well the detector can distinguish events with single elec- 
tron induced and two photon induced electromagnetic 
(e.m.) showers. Backgrounds due to — > v T conver- 
sion, v T charged current reactions in the detector, and 
subsequent t~ — > e~v e u T decays are small because of 
the low energy of the beam and the consequent low cross 
section for r production [14]. 

Because of the rapid fall in the BVLB neutrino spec- 
trum beyond 4 GeV (Figure 3), the largest contribu- 



tion to the 7T° background will come from neutral cur- 
rent events with a single ir° in the final state. It is well 
known that resonant single pion production in neutrino 
reactions has a rapidly falling cross section as a function 
of momentum transfer, q 2 , [26] up to the kinematically 
allowed value. This characteristic alone suppresses this 
background by more than 2 orders of magnitude for ir° 
(or shower) energies above 2 GeV [14]. Therefore a mod- 
est 7T° background suppression (by a factor of ~ 15 below 
2 GeV and ^2 above 2 GeV) is sufficient to reduce the 
7T° background to manageable level over the entire spec- 
trum. The electron neutrino contamination in the beam 
from decays of muons and kaons, is well understood to be 
0.7% of the muon neutrino flux [27], with a similar spec- 
trum (Figure 3). The experimentally observed electron 
neutrino energy spectrum will therefore have three com- 
ponents: the rapidly falling shape of the tt° background, 
the spectral shape of the v e beam contamination slightly 
modified by oscillation, and the oscillatory shape of the 
appearance signal. The shape of the appearance spec- 
trum will be well known because of the precise knowledge 
of Am 2 2 from the disappearance measurement as well as 
the improved knowledge of Am 2 i from KamLAND [5]. 
These distinguishing spectra will allow experimental de- 
tection of — > v e with good confidence. Figure 6 shows 
a simulation of the expected spectrum of reconstructed 
electron neutrino energy after 5 x 10 7 sec of running. The 
parameters assumed are listed in the figure as well as the 
caption. 
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FIG. 5: (color) Resolution including statistical and system- 
atic effects at 90% confidence level on Am 2 ^ and sin 2 2023 for 
the 2540 baseline experiment; assuming 1 MW, 0.5 MT, and 
5 x 10 7 sec of exposure. We have included a 5% bin-to-bin 
systematic uncertainty in the energy calibration as well as a 
5% systematic uncertainty in the normalization. We have not 
included a systematic uncertainty on the global energy scale; 
this should be added in quadrature to the expected resolution 
on Am.32. The expected resolution from the MINOS experi- 
ment at Fermilab (dashed) and the allowed region from Super 
Kamiokande (dotted) is also indicated. 



Figure 6 further illustrates the previously described 
three regions of the appearance spectrum: the high en- 
ergy region (> 3 GeV) with a matter enhanced (for NO) 
appearance has the main contribution to the background 
from v e contamination in the beam; the intermediate re- 
gion (1 — > 3 GeV) with high sensitivity to the CP phase, 
but little dependance on mass ordering, has approxi- 
mately equal contribution from both background sources 
(y e and 7r°); and the low energy region (< 1.5 GeV) , 
where the effects of the CP phase and Amfi dominate, 
has the main background from the 7r° events. Matter- 
enhancement of the oscillations has been postulated for 
a long time without experimental confirmation [12]. De- 
tection of such an effect by observing a matter enhanced 
peak around 3 GeV will be important. However, in the 
case of RO mass ordering this enhanced peak will be 
missing, but the effect (depending on 8cp) on the rest of 
the electron neutrino spectrum will be small. 



The very long baseline combined with the low energy 
spectrum make it posible to observe — > v e conversion 
even if sin 2 2#i3 = because of the contribution from 
Am^ if the solar neutrino large mixing angle solution 
(LMA) holds. In Figure 7 we show the expected appear- 
ance spectrum for Am^ = 7.3 x \0~ b eV 2 (LMA-I) with 
90 excess events and Am| x = 1.8 x lQ- 4 eV 2 (LMA-H) 
with 530 excess events above background. The LMA- 
I and LMA-II solutions with their respective confidence 
levels are discussed in [28] . This measurement will be sen- 
sitive to the magnitude and knowledge of the background 
because there will be no oscillating signature to distin- 
guish the signal. We estimate that the statistical and 
systematic errors on this measurement will allow us to 
determine the combination Amjj x sinWyi to a precision 
of about 10% at the LMA best fit point. This is compet- 
itive with the expected measurement from KamLAND. 
However, it will be done in the explicit (y^ — > v e ) ap- 
pearance mode, which is qualitatively different from the 
measurements made so far in the SNO or the KamLAND 
experiments. Therefore, such a measurement will have 
important implications on the search for new physics such 
as sterile neutrinos. 

In Figure 8 we show the 90 and 99.73% (3 sigma) 
C.L. sensitivity of the BVLB experiment in the variables 
sin 2 2^i3 versus Sep- The actual limit obtained in the 
case of a lack of signal will depend on various ambigui- 
ties. Here we show the 99.73% C.L. lines for NO and RO, 
on the right hand of which the experiment will observe an 
electron appearance signal due to the presence of terms 
involving #13 with greater than 3 sigma significance and 
thus determine the corresponding mass ordering. Cur- 
rently from atmospheric data there is a sign uncertainty 
on 623 = ±7r/4; this introduces an additional ambiguity 
onto Figure 8 of Sep — > 5c p + 7T. For this plot we have 
assumed that the other parameters are well known cither 
through other experiment or by the disappearance mea- 
surement: Am^ = 7.3 x 10~ 5 eV 2 , Ato 2 9 = 0.00 25 eV 2 . 
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The values of sin 2 20 12 and sin 2 26> 2 3 are set to 0.86, 
1.0, respectively. The sensitivity to sin 2 29i 3 is approx- 
imately constant in the Am| 2 range allowed by Super 
Kamiokande [14]. The value of Amfr affects the modu- 
lation of the sin 2 26*13 sensitivity with respect to Sep- If 
there is no excess electron appearance signal, other than 
the expected signal due to Am^ alone, then it would 
indicate a switch to anti-ncutrino running to determine 
if the RO hypothesis with parameters on the left hand 
side of the dashed line in Figure 8 applies [20] . 
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FIG. 6: (color) The expected electron neutrino spectrum for 
3 values of the CP parameter Sep including background con- 
tamination. The points with error bars are for 5c p = 135°; 
the error bars indicate the expected statistical error on each 
bin. The histogram directly below the error bars is for Sep = 
45° and the third histogram is for 8c p = —45°. The hatched 
histogram shows the total background. The v e beam back- 
ground is also shown. The plot is for Am§ 2 = 0.0025 eV 2 . We 
have assumed sin 2 26>i 3 = 0.04 and Am|i = 7.3 x 10~ 5 eV 2 . 
The values of sin 2 2(9i2 and sin 2 2#23 are set to 0.86, 1.0, re- 
spectively. Running conditions as in Figure 4. 



Sensitivity to the CP Violation Parameter 

To get a qualitative understanding for the measure- 
ment of CP violating parameters we compare the size 
of terms involving Sep with the first term in Equation 
1. This ratio is seen to be approximately proportional 
to asmScpi^m^L/AE^). As shown in Figure 2, the 
fractional contribution from the CP violating terms in- 
creases for lower energies at a given distance. The energy 
dependence of the CP effect and the matter effect tend 
to be opposite and therefore can be distinguished from 
each other from the energy distribution using neutrino 
data alone. On the other hand, the statistics for a given 
size detector at a given energy are poorer by one over 
the square of the distance, but the term linear in sin<5cp 
grows linearly in distance [13]. The statistical sensitiv- 
ity (approximately proportional to the square root of the 
event rate) to the effects of CP violation, therefore, is 
independent of distance because the loss of event rate 



FIG. 7: (color) The expected spectrum of electron neutrino 
events for sin 2 29is — 0. Important parameters in this figure 
are Am^ = 7.3 x 10" 5 eV 2 (LMA-I) or Am^ = 1.8 x 10~ 4 
(LMA-II) and sin 2 2#i2 = 0.86. All other parameters and the 
running condition as in Figure 4. 



and the increase of the CP effect approximately cancel 
each other in the statistical merit if backgrounds remain 
the same. Therefore, the two important advantages of 
the BVLB approach are that the CP effect can be de- 
tected without running in the anti-neutrino mode and 
the sensitivity to systematic errors on the background 
and the normalization is considerably reduced because 
the fractional size of the CP effect is large. For exam- 
ple, in Figure 2, the CP effect at Sep = 7t/4 in the first 
oscillation peak is ~ 20% while the effect in the second 
oscillation peak is more than 50%. Therefore, it is unnec- 
essary to know the background and the normalization to 
better than 10% to obtain a significant measurement of 
Sep at the second oscillation maximum. There could be 
a contribution to the systematic error from the theoreti- 
cal calculation of the probability shown in Figure 2. Be- 
cause of the very long baseline, this probability depends 
on the Earth's density profile which is known to about 
5%. Random density fluctations on that order will lead 
to a relative systematic uncertainty in the v e appearance 
probability of about 1% [29], which is not significant for 
the BVLB method, but could be significant in the case 
of an experiment that performs the CP measurement at 
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FIG. 8: (color) 90 and 99.73% CL. contours of the BVLB 
experiment in the variables sin 2 2#i3 versus 5c p for the natu- 
ral (NO) and reversed ordering of parameters (RO). The solid 
lines are for NO, the left line for 90% and right line for 99.73% 
CL. The dashed and dotted lines are for 90 and 99.73% CL. 
for RO. On the right hand side of the 99.73% CL. lines the 
experiment will observe an electron appearance signal with 
greater than 3 sigma significance. 



not change rapidly with respect to 5c p, therefore the 
resolution on Sep becomes poor. 

Because of the multiple node spectrum the appearance 
probability is sensitive to both the sin 5c p and cos 5c p 
terms in Equation 1. This eliminates the 5c p — > (it — 
6cp) ambiguity normally present in a narrow band single 
node experiment [8]. The sign uncertainty #23 = ±7r/4 
introduces a 5c p — ► (5c p + i") ambiguity. If Am^ is 
known poorly, there will be a contribution to the 5c p 
resolution due to the correlation between 5c p and Ato 2 ^. 
We expect to know Ato^i from KamLAND with precision 
of < 10% [5] ; this does not introduce a dominant source of 
uncertainty on 5cp in our experimental method. If a four 
parameter, Am^, sin 2 2#i2, sin 2 26*13, and 5c p, fit to the 
data is performed, the correlations between the param- 
eters dilute the sensitivity to each parameter. However, 
the ambiguities and correlations do not reduce the ability 
of the experiment to establish that 3-generation neutrino 
mixing contains a non-zero complex phase and hence a 
CP violating term. This is seen if we consider the resolu- 
tion on the quantity Arn^j x Jcp which is CP violating 
and by definition does not exhibit any of the above cor- 
relations. The measurement of this quantity will simply 
depend on the statistics in the low and medium energy 
components of the spectrum. 



the first oscillation maximum with a shorter baseline. 

An experiment with a wide band beam is needed to ex- 
ploit the energy dependence of the CP effect to extract 
the value of 5c p and sin 2 20 13 using only the neutrino 
data. The energy dependence also reduces the correla- 
tion between the two parameters. Figure 9 shows the 
expected resolution on sin 2 2#i3 versus 5c p at a particu- 
lar case, sin 2 26*13 = 0.04 and 5cp = 7r/4 with all other 
parameters fixed as indicated. The la error on 5c p from 
this measurement as a function of Sep is shown in Figure 
10. It is clear that the precision of the measurement will 
be limited if sin 2 26*13 < 0.01 because of the background. 
However, for larger values the resolution on 5c p is ap- 
proximately independent of sin 2 2#i3 . This result can be 
understood by examining Equation 1 in which the first 
term grows as sin 2 2#i3 and the CP violating term grows 
as sin 26*13. Therefore the statistical sensitivity to the CP 
violating term, which is the ratio of the second term to 
the square-root of the first term, is indenpendent of 6*13 
as long as the background docs not dominate [13]. At 
Sep ~ 135° the event rate reaches maximum and does 



Comparison to Other Techniques 

The alternatives to the BVLB plan fall into two cat- 
egories: experiments with conventional horn focused 
beams over baselines of several hundred km and with 
beams based on neutrinos from muon decays in a high 
intensity muon storage ring. For a recent review of muon 
storage ring based neutrino factory physics as well as 
comparisons to conventional neutrinos beams see [30]. 

Recently, for the conventional beam an "off-axis" beam 
is preferred because the energy spectrum is narrower and 
can be tuned, within a limited range from 1 to 2 GeV by 
choosing the off-axis angle [6] . A narrow energy spectrum 
could be better suited for identifying — > v e conversion 
because the neutral current backgrounds from interac- 
tions of high energy neutrinos are lowered. The JPARC 
to Super Kamiokandc project (baseline of 295 km) has 
adapted this strategy and a new initiative at Fermilab 
based on the NuMI bcamlinc (baseline of 732 km) is also 
considering it [7]. Both of these projects focus on the ob- 
servation of a non-zero value for #13 . Since it is not pos- 
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Am§ 2 = 0.0025 eV 2 , and Amji = 7.3 x 10~ 5 eV 2 . The values 
of sin 2 2#i2 and sin 2 2023 are set to 0.86, 1.0, respectively. 
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FIG. 10: (color) 1 sigma error on 8c p in degrees for sev- 
eral values of sin 2 2#i3 . All other parameters are held fixed 
assuming they will be known before this measurement. 



sible to eliminate the ^0.5% v e contamination in a con- 
ventional beam (even if it is "off-axis" ) the sensitivity 
of all projects based on a conventional beam, including 
the BVLB plan, is limited to sin 2 2#i3 ~ 0.005 (ignoring 
the effect of Sep)- Sensitivity below this value is back- 
ground limited and requires very good understanding of 
the systematic errors associated with the purity of the 
beam. The BVLB plan does, however, extend this sensi- 
tivity to values of Am 2 down to 0.001 eV 2 . The "off-axis" 
projects must be upgraded to several Megawatts of beam 
power and detectors of several hundred kT (the second 
phase of the JPARC project calls for 4 MW beam and a 1 
megaton detector) to perform the CP measurement. The 
narrow band nature of the beam requires running in the 
anti-neutrino mode to observe CP violation. The BVLB 
plan uses a wide band beam over a much longer distance 
to achieve the same sensitivity to sin 2 2#i3 over a wider 
range of Am| 2 . It also has the ability to see a multiple 
node structure in both the disappearance and appear- 
ance channel; thus making the observation of — > v e 
experimentally more robust against unanticipated back- 
grounds. In addition, the unique aspect of the wide band 
beam over the very long baseline is the ability to measure 
Scp with few ambiguities using neutrino data alone. 

The muon storage ring based "neutrino factory" solves 
the background problem faced with conventional beams. 
A and a v e is produced in each (i + decay. An ap- 
pearance search looks for the conversion of v e into by 
looking for events with negatively charged fi~ particles in 
the presence of charged current events with positively 
charged \i + particles in the final state. Such a physics 
signature should potentially have very low background. 
Therefore, a neutrino factory based v e — » search is 
estimated to be sensitive to sin 2 26*i3 perhaps as small 

as 10~ 4 10~ 5 . Such a facility is very well tuned for 

measurement of sin 2 26>i3. Nevertheless, because of the 
relatively narrow energy spectrum and higher energy of 
the neutrino factory, it lacks the ability to observe mul- 
tiple nodes in the observed spectrum. For the neutrino 
factory the effects of CP violation are observed by com- 



paring v e 



and 



in two different runs of ji 



and n~ decays in the presence of very large matter ef- 
fects. The narrow spectrum also limits the ability of the 
neutrino factory to resolve the ambiguities in the 5c p 
determination regardless of the value of sin 2 26*13. This 
has been extensively discussed in the literature [8, 25] 
and one of the proposed solutions is to construct sev- 
eral detectors with different baselines to uniquely deter- 
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mine 5c p versus sin 2 26*13. Although the neutrino factory 
can extend the reach for sin 2 26*13 to very small values 
because the backgrounds to the appearance search are 
small, the resolution on 5cp - disregarding the ambigu- 
ities inherent in the neutrino factory approach - is com- 
parable to the BVLB resolution. We therefore conclude, 
that if sin 2 29\3 > 0.01 the BVLB experiment is either 
equal to or exceeds the physics reach for neutrino mixing 
parameters of the neutrino factory. 

Other suggestions in the literature include v e or v e 
beams from either decays of radioactive nuclei acceler- 
ated to high energies (beta-beam) [31] or P e 's from nu- 
clear power generation reactors such as the beam used 
for KamLAND [5]. Both of these concepts aim to cre- 
ate very pure, but low energy beams. For CP sensitivity 
both v e and v e must be produced in separate runs in the 
beta-beam concept. The low background should enable 
the beta-beam concept to have good sensitivity to #13 , if 
sufficient flux can be obtained. The reactor beams are be- 
low charged current reaction threshold and therefore 
must rely on disappearance to search for a non-zero value 
of #13. A disappearance based search for 6*13 will most 
likely be limited by systematic errors. Further comments 
on these techniques must await more detailed studies. 

Conclusions 

We have simulated and analyzed a feasible, very long 
baseline neutrino oscillation experiment consisting of a 
low energy, wide band neutrino beam produced by 1 MW 
of 28 GeV protons incident on a carbon target with mag- 
netic horn focusing of pions and kaons, and a 500 kT 
detector at a distance of > 2500 km from the neutrino 
source. The BVLB neutrino beam with a total intensity 
of about 4.7 x 10~ b v/m 2 /POT at a distance of 1 km 
from the target could be provided by an upgrade to the 
BNL-AGS [17]. 

The single BVLB experiment could produce measure- 
ments of all parameters in the neutrino mixing ma- 
trix through observation of the disappearance channel, 
v n ~ * v ii an d t ne appearance channel, — > v e . The 
experiment is also sensitive to the mass ordering of neu- 
trinos using the observation of the matter effect in the 
appearance channel through the currently unknown pa- 
rameter sin 2 26*13. Nevertheless, the experiment is in- 
tended primarily to measure the strength of CP invari- 
ance violation in the neutrino sector and will provide a 



measurement of the CP phase, 5c p or alternatively the 
CP violating quantity, Jcp , if the one currently unknown 
neutrino oscillation mixing parameter sin 2 26 \$ > 0.01, a 
value about 15 times lower than the present experimen- 
tal upper limit. We point out that for a given resolu- 
tion on 5c p the number of neutrino events needed which 
determines the detector size and beam intensity, is ap- 
proximately independent of the baseline as well as the 
value of sin 2 20 13 as long as the electron signal is not 
background dominated. Therefore, the concept of very 
long baseline (> 2000 km) is attractive because it pro- 
vides access to much richer physics phenomena as well 
as reduces the need to understand systematic errors on 
the flux normalization and background determination. A 
shorter baseline (< 1500 km) will obviously limit the 
reach of the experiment for Arn^i x sin26>i2 in the ap- 
pearance mode as well as reduce the resolution on the CP 
parameter. On the other hand, a much longer baseline 
(> 4000 km) will result in a matter effect that is large 
enough to dominate the spectrum and make the extrac- 
tion of the CP effect more difficult. The larger distance 
will also make it necessary to make the neutrino beam 
directed at a higher angle into the Earth, a technical chal- 
lenge that may not be necessary given the current values 
of neutrino parameters. Lastly, we have also shown that 
most of this rich physics program including the search for 
CP effects can be carried out by neutrino running alone. 
Once CP effects are established with neutrino data, anti- 
neutrino data could be obtained for more precision on 
the parameters or search for new physics. The shape of 
the disappearance signal over multiple oscillations from 
neutrino running alone (as well as in combination with 
anti-neutrino running) can be used to contrain effects of 
new physics from: sterile neutrino mixing, extra dimen- 
sions, exotic interactions with matter, etc. 



It has not escaped our notice that the large detector 
size (>500 kT), mandated by the above described neu- 
trino program, naturally lends itself to the important 
physics of proton decay, supernova detection, etc. Al- 
though some experimental requirements such as shielding 
from cosmic rays, detector shape, or photo-sensor cover- 
age, may differ, it is clear that they can be resolved so 
that this large detector can also be employed to address 
these frontier physics issues. 
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